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ABSTRACT 



It is proposed that the observed near-surface inflows towards the active regions and sunspot zones provide a nonlinear feedback 
mechanism that limits the amplitude of a Babcock-Leighton-type solar dynamo and determines the variation of the cycle strength. 
This hypothesis is tested with surface flux transport simulations including converging latitudinal flows that depend on the surface 
distribution of magnetic flux. The inflows modulate the build-up of polar fields (represented by the axial dipole) by reducing the tilt 
angles of bipolar magnetic regions and by affecting the cross-equator transport of leading-polarity magnetic flux. With flux input 
derived from the observed record of sunspot groups, the simulations cover the period between 1874 and 1980 (corresponding to 
solar cycles 1 1 to 20). The inclusion of the inflows leads to a strong correlation of the simulated axial dipole strength during activity 
minimum with the observed amplitude of the subsequent cycle. This in agreement with empirical correlations and in line with what 
is expected from a Babcock-Leighton-type dynamo. The results provide evidence that the latitudinal inflows are a key ingredient in 
determining the amplitude of solar cycles. 
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1. Introduction 

The correlation between the amplitude of the solar polar field 
around a ctivity minimum and the strength of th e subsequent cy- 



cle (e .g., ISchatten et al. 



^^^^^ „ , - JigTSHCh oudhuri. 20081: ICameron et al.L 

see, however. iLay den et all |1991|) suggests that the Sun's 
axial dipole field as observed around activity minimum re- 
flects the source for the generation of the toroidal magnetic flux 
of the subsequent cycle, as opposed to being a mere epiphe- 
nom enon of an otherwise fully hidden subsurface dynamo pro- 
cess (ICameron & Schussleill2007l : ISc huss"ler. 2007). Such a role 
of the polar field is a key feature of the Ba bcock-Leighton (BL) 
dynamo model (see ICharbonneaui 1201 Oi for a review of so- 
lar dynamo theory). In this type of dynamo, the poloidal mag- 
netic field results from the systematic tilt (with respect to the 
azimuthal/longitudinal direction) of sunspot groups and bipolar 
magnetic regions in combination with the subsequent redistribu- 
tion of their magnetic fluxes by near-surface flows. The crucial 
mechanism for the reversal of the Sun's global dipole field is 
the preferred transport of preceding -polarity flux of bipolar mag- 
netic regions across the equ ator into the other solar hemisphere 
(ICameron & Schussleiil20(j7,) . which leads to the reversal of the 
polar fields and the build-up of an axial dipole field of opposite 
polarity. Toroidal flux is then generated by the action of differ- 
ential rotation winding up poloidal field lines threading the Sun. 

So far we have described the BL dynamo as a linear pro- 
cess. In addition, at least one nonlinear feedback mechanism is 
required that limits the amplitude of the generated magnetic field 
and controls the variation of the cycle amplitude. In the frame- 
work of the BL dynamo, potential feedback mechanisms include 
the back-reaction of the magnetic field on 1) differential rotation, 
2) the tilt angles of sunspot groups, and 3) cross-equator trans- 
port of magnetic flux. 



Concerning possibility 1), the only empirically well- 
established cycle-related variation of the differential rotation 
are zonal flows of a few m s which at lo w latitud es ar e 
associated with the activi ty belts dHovyard & Labonte. '198' 
IZhao & Kosovichevll2004tlGizonll2004l: lGizon & Birch, 2001 
Clear evidence for a more global reaction on differential rotation 
that would affect the generation of the toroidal magnetic field 
has not been obtained so far. 

As to mechanism 2), the back-reaction on the tilt angles of 
sunspot groups, Dasi-Espuig et al. (2010) found cycle-to-cycle 
variations of the average tilt angle, which are anti-correlated with 
the cycle strength: strong cycles show smaller average tilt an- 
gles than weak cycles. Since the tilt angles affect the strength 
of the source for the poloidal field in the BL dynamo, their an- 
ticorrelation with cycle strength provides a p ossible nonlinear 
saturation mechanism. iCameron et alJ (120101) included the ob- 
served tilt-angle variations in a surface flux transport (SFT) sim- 
ulation with sources taken from the combined Royal Greenwich 
Observatory (RGO) and the Solar Optical Observing Network 
(SOON) sunspot records. Considering the period for which tilt- 
angle data are available (1913-1986), they found that the max- 
imum of the polar field (occurring around activity minimum 
of a given cycle) in the SFT simulation is correlated with the 
observed strength of the subsequent activity cycle. This sug- 
gests that the anti-correlation of the average tilt angle with cycle 
strength contributes to the variations of the cycle amplitudes and 
thus to the nonlinear feedback by which the solar dynamo satu- 
rates at the observed levels. 

Which physical process could lead to the observed anti- 
coiTelation? One possibility is a change in the strength or ef- 
fectiveness of the Coriolis force that tilts magnetic flux tubes 
rising in the convection zone (e.gjDlSilva & Choudhuri, 11991 
iFan et all Il994t ICafigari et all Il995[ll998i iFani i2009) . In the 
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absence of sufficient observational information about tlie sub- 
surface structure and dynamics of the magnetic field, it is 
difficult to evaluate this mechanism. Another possible pro- 
cess is more amenable to observation and quantitative test: 
tilt-angle changes caused by the redistribution of magnetic 
flux after emergence by near-surface flows resulting from 
the presence of the magnetic field. A well-observed exam- 
ple of such flows are the inflows towards act ive regions and 
the activity belts (e.g., iHaber et all |2 002; Zha o & Kosovichevl 
200l'Gizon'.'200?.'Gizon & Re mpelLl200&,Svanda et aLl lIOOSt 
Gon zalez Hernandez et al.. 201)3), possibly resulting from a 
temperature deficit due to the exces s radiance of small-scale 
magn etic flux concentrations jSpruitl 120031: iGizon & Rempell 
l2008h . Collectively, these inflows drive mean flows towards the 
activity belts and zonal flows, whose amplitudes depend on the 
amount of magnetic flux present in active regions, i.e., on the 
overall strength of a cycle. The transport of active-region flux by 
both the converging latitudinal flows and the longitudinal zonal 
flows tends to reduce the tilt angles of the coiTesponding bipo- 
lar magnetic regions. This mechanism was recently studied by 
[jiang et al. (2010) with SFT simulations including an ad-hoc 
model of the inflows in terms of latitudinal flows towards the 
activity belts. Considering individual bipolar regions as well as 
synthetic solar cycles, they found that the inflows reduce the lat- 
eral separation of the magnetic polarities, leading to a weaken- 
ing of the contribution of a given bipolar regio n to the polar field. 
Building on these results, iJiang et al] (1201 lalfB) included cycle- 
dependent tilt angles in their reconstruction of the evolution of 
the surface flux and open heliospheric flux since 1700. 

There is another aspect of the converging inflows towards 
active regions, which is related to possibility 3) for a nonlin- 
ear effect listed further above: a feedback on the cross-equator 
transport of magnetic flux. For active regions not too far from 
the equator, the inflows often extend over the equator (see, e.g.. 
Fig. 18 in Gizon et all 120101) . which enhances the otherwise 
purely diffusive cross-equator transport of preceding-polarity 
flux and, therefore, tends to increase the resulting polar field. 
This effect is weaker for strong cycles because then active re- 
gions start to emerge at higher latitudes and the activity maxi- 
mum is reached earlier (the Waldmeier effect), at a time when 
the activity belts are farther away from the equator. Therefore, 
this represents a negative feedback on the growth of the polar 
field. 

In this paper, we present results of SFT simulations includ- 
ing mean latitudinal flows towards the activity belts, which were 
determined self-consistently from the simulated distribution of 
magnetic flux, thus introducing a nonlinear element into the SFT 
model. Using the RGO/SOON sunspot group data to determine 
the flux in emerging bipolar magnetic regions, we followed the 
evolution of the surface flux from 1874 until 1980. In order to 
evaluate the possible relevance of the nonlinearity for determin- 
ing the cycle strength in the framework of a BL dynamo, we 
considered the correlation between the maximum of the mag- 
netic flux associated with the axial dipole component in a given 
cycle resulting from the SFT simulation with the observed am- 
plitude of the subsequent cycle. 

The advantage of the approach chosen here is that it deals 
exclusively with observationally constrained near-surface pro- 
cesses and thus does not require a complete dynamo model in 
order to evaluate the potential of the nonlinearity. While the sur- 
face evolution addressed by SFT simulations is constrained by 
observations, a full dynamo model involves assumptions of the 
subsurface dynamics, which at present are essentially uncon- 
strained. 



2. Surface flux transport simulation 

2.1. Transport equation 

We used the SFT code of Baum ann et al.l (|2004|) in the version 
described bv^ Jiang et al.i (,20 10) to simulate the evolution of the 
magnetic flux at the solar surface under the influence of large- 
scale surface flows (differential rotation and meridional flow) 
and a diffusivity representing the random motion of the magnetic 
flux elements due to the (numerically unresolved) supergranular 
flows. The magnetic field at the solar surface is assumed to be 
radial. The governing equation of the SFT model is 

^ = - oj{A)^- ^ ^ [\v{A) + ^v{AJ)\Bco^A] 
dt d(p RQCOsAdAv^ J I 

I d I dB\ 1 d^B 

+ VH T— COS/l— H — r- 

'"[RQCQ?.AdA\ dAj Rlco&^Adf 
+ SiA,ip,t) + D(7j,), (1) 

where B{A, (f>) denotes the radial component of the magnetic 
field, A is the heliographic latitude, and (p is the heliographic 
longit ude. The express ions for the synodic differential rotation 
(from lSnodgrassl Il983l) . 

Lo(A) = 13.38 - 2.30 sin^^- 1.62 sin"* /I [7day], (2) 

and for the time-independent part of the meridional flow (cf. 
van Ballegooiien et al., 1998), 

y(^) _ f 1 1 sin(;ri/75°) m s"' where M |< 75° 
\0 otherwise, 

as well as the choice of the horizontal and radial diffusivities 
(lBaumanneta l.','2006') , tjh = 250 km ^s'^ and D = 0, were the 
same as used in Camer on et alJ (|2010|). The velocity Au(A, t) in 
Eq. ([TJ represents the longitude-averaged inflows toward the ac- 
tivity belts, which are described by a magnetic-field-dependent 
local perturbation of the meridional flow; its construction is ex- 
plained in Sec. 12.31 This term renders the model nonlinear. 

2.2. Flux emergence 

The source term, S {A, (p, f), in Eq. ([TJ represents the emergence 
of magnetic flux at the surface in the form of bipolar mag- 
netic regions appearing at the latitudes and longitudes of sunspot 
groups taken from the combined RGO/SOON sunspot record 
dBalmaceda et ill l2009h at the time of their maximum surface 
area. The sunsp ot group size wa s conve rted to magnetic flux flux 
as described in ICameron et al. and references therein), 

the calibration being based upon the total unsigned magnetic 
flux derived from the synoptic magnetic map^ taken with the 
Michelson Doppler Imager (MDI) instrument onboard of the 
Solar Heliospheric Observatory (SoHO) spacecraft. We assumed 
that the latitude dependence of the initial tilt angle, a, of the 
new bipolar magnetic regions has the form a - T yfJ, using 
T = 1.42 for all cycles (angles measured in degrees), so that 
cycle-dependent changes of the tilt angles result solely from the 
magnetic-field dependent inflows toward the activity belts. 

The transition from the RGO data covering the period up 
to 1976 to the subsequent SOON data involves cross-calibration 
issues. For instance, the minimum area required for a sunspot 
group to be included changed from 1 millionth of a hemi- 
sphere in the RGO data to 10 millionths of a hemisphere in the 

' ihttp : //soi ■ Stanford . edu/magnetic/index6 . html | 
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Fig. 1. Time-latitude diagrams of various quantities from a SFT simulation run with magnetic-field dependent meridional inflows to- 
wards the activity belts. From top to bottom: longitudinally averaged signed radial magnetic field, longitudinally averaged unsigned 
radial magnetic field, meridional flow velocity, and meridional flow perturbation (Av{A, t), saturated for values exceeding ±2 m-s '). 
Black vertical lines indicate the termination of the RGO sunspot data in 1976. 



SOON data. In order to account for this change. lBalmaceda et alj SOON data to determine a correction factor of approximately 
(.2009.) used other datasets overlapping both the RGO and 1.4 for the sunspot group ai^eas in the SOON data (see also 
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IWilson & Hathawavl l2005h . While a global correction factor 
might be adequate for the total sunspot areas, applying it uni- 
formly to individual sunspot groups of all sizes (as used in the 
source term of SFT simulations) is hardly justified since it prob- 
ably grossly overestimates the areas of the large sunspot groups. 
This affects particularly strongly the variation of the polar field, 
to which the large sunspot groups contribute most because of 
their large latitude separations and fluxes. Since the overesti- 
mated sources lead to unrealistic results for the polar fields (and 
the closely related axial dipole component), which are at our fo- 
cus of interest in this paper, we consider the SFT results only un- 
til the activity minimum between cycles 20 and 21, which nearly 
coincides with the termination of the RGO data in 1976. 

2.3. Latitudinal infiows 

We treated the cumulative effect of the inflows toward active re- 
gions in terms of a local perturbation, Au{A, t), of the axisym- 
metric meridional flow. The amplitude of the background merid- 
io nal flow wa s kept constant. This is consistent with the result 
of ICameron & Schiissler (2010), who showed that the active- 
region inflows fully explain the variations in the low-degree 
spher ical -harmonic decomposition of the observed meridional 
flow (iHathawav & Rightmire, 2010), without requiring a change 
in the overall meridional flow amplitude. A time-independent 
gl obal meridional ci rculation is also consistent with the finding 
bv lHathawavl (|201 ih that the latitudinal drift rates of the activity 
belts are the same for all cycles when measured with respect to 
the epochs of the cycle minima. The perturbation of the merid- 
ional flow in our SFT simulation s was therefore modeled usin g 
almost the same procedure as in ICameron & Schiissled (l2010h . 
Taking the instantaneous latitude profile of the azimuthally av- 
eraged unsigned radial magnetic surface field, 

{\B\){A,t)^-^ r |B(0,^,O)|d0, (4) 
Jo 

we set the inflow speed towards the activity belts proportional to 
the smoothed derivative of {\B\) with respect to latitude, viz. 

where the choice of cr effects a smoothing in latitude with 
a full width at half maximum of 20°. The prefactor cp - 
9.2 ms 'G 'deg was calibrated in ICameron feSchiissleil (120101) 
by requiring that the amplitude of t he infl ow should be compara- 
ble to that reported by Gizon et aTl (|201tf) from helioseismic ob- 
se rvations during cycle 23. The only difference to the procedure 
of ICameron & Schiissled (1201 Ol) is the factor (cos A/ cos 30°) in 
Eq. (jSj, which suppresses unrealistically strong flow perturba- 
tions that would otherwise result at high latitudes from the gra- 
dient of the polar fields. 

For a nonlinear SFT simulation including the effect of 
magnetic-field-dependent inflows, Fig.[T]shows time-latitude di- 
agrams of the longitude-averaged radial magnetic field (signed 
and unsigned) as well as of the meridional flow speed and its per- 
turbations, Au, due to the modeled inflows. Stronger cycles are 
associated with bigger perturbations of the meridional flow. For 
very high cycle amplitude, the flow direction may even locally 
reverse, which can be seen most prominently during cycle 19. 
Depending on the asymmetry between the magnetic field distri- 
butions in the two hemispheres, cross-equator flows occasionally 
arise. 
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Fig. 2. Magnetic flux corresponding to the axial dipole, <I> dipole, 
from the SFT simulation with inflows (reference case, black 
curve), without inflows (red curve), and with strong inflows (am- 
plitude multiplied by 1.5 with respect to the reference case, blue 
curve). The black vertical line indicates the termination of the 
RGO sunspot data in 1976. 

3. Results 

3.1. Effects of the inflows on the axial dipole 

The latitudinal inflows towards the activity belts affect the near- 
surface evolution of the magnetic field in a nonlinear fashion, 
since they depend on the magnetic field distribution and strength. 
To investigate the potential of this effect as a nonlinear feedback 
in flux-transport dynamos, we compare SFT simulations with 
and without inflows. A relevant quantity for the amplitude of 
flux-transport dynamos is the maximum axial dipole component 
of the surface field around activity minimum. This field repre- 
sents the main part of the poloidal field from which the toroidal 
field of the next activity cycle is generated by differential rota- 
tion. Note that the axial dipole component is axisymmetric and 
therefore independent of the presence of zonal flows since any 
differential rotation does not affect the longitude-averaged rad ial 
surface field in a SFT simulation (e.g.. lBaumann et al.Ll2004l) . 

The time development of the magnetic flux associated with 
the axial dipole, Odipoie, for cycles 12 through 20 is shown in 
Fig.|2] Three SFT simulation runs are compared: the 'reference' 
case (cf. Fig.[T]) with calibrated inflows towards the activity belts 
(black line), the case without inflows (red line), and a case for 
which the inflows were enhanced by a factor 1.5 with respect to 
the reference case (blue line). The first cycle in the simulation 
is not considered further sin ce it is still dominated by the ini- 
tial condition (for details, see lCameron et alll2010t) . In both the 
reference case and the case with no inflows, Odipoie peaks at or 
shortly before the activity minima. 

For all cycles except cycles 19, we see that the inflows in the 
reference case tend to increase Odipoie- We can understand this 
result by considering that the preferential cross-equator trans- 
port of leading-polarity magnetic flux in tilted bipolar mag- 
netic regions is the essential effect that produces a flux im- 
balance in each hemisphere, leading to the polarity reversal 
and build-up of the o pposite-polarity axial dipole field (e.g., 
ICameron & Schiissler. 2007). There are two opposing effects of 
the latitudinal inflows on the cross-equator transport: 1) for low- 
latitude active regions, the inflows are broad enough to cross 
the equator and thus enhance the flux transport from the other 
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curve) leads to a weakened axial dipole after cycle 19 and a much 
stronger dipole after cycle 20. 

We studied the sensitivity of the axial dipole on the strength 
of the inflows by running SFT simulations with inflow ampli- 
tudes multiplied by factors of 0.5 and 1.5, respectively, with re- 
spect to the reference case. While in the case of weakened in- 
flows the time evolution of Odipoie is essentially the same as in 
the case of no inflows (red curve in Fig. |2]i, for the enhanced 
inflows the deviations from the reference case are drastic (blue 
curve in Fig.|2|. The overly strong reduction of the cross-equator 
flux transport in this case leads to some almost failed polar field 
reversals and a resulting strong 22-year periodicity of <!' dipole- 

For a comparison of the SFT models with empirical results, 
Fig.[3]shows the relation between the maxima of Odipoie from our 
SFT simulations and the open heliospheric magnetic flux dur- 
ing these perio ds determined from the v ariations of the geomag- 
netic flfl-index dLockwood et al.L [2009 '). We may safely assume 
that the open magnetic flux during activity m inimum periods is 
stron gly dominated by the axial dipole (e.g.. lWang & Sheelevl 
12009 ). The upper panel of Fig. [3] gives the result for the sim- 
ulation without inflows, which shows a weak (r = 0.55) and 
insignificant {p = 0.15) correlation; there is also a clear devi- 
ation of the regression line (solid) from the bisectrix (dashed). 
In contrast, for the reference case with inflows (lower panel of 
Fig.O we find a good correlation (r = Q.lA,p - 0.022) and a 
regression line nearly coinciding with the bisectrix. It is obvious 
from the blue curve in Fig.|2]that the case with enhanced inflows 
is inconsistent with the empirically determined open fluxes. 




n L ^ , , I , , , I , , , I , , , 

2 4 6 8 

*aipo. ["^O" Wb] 




Fig. 3. Maxima of Odipoie from SFT simulations vs. open flux 
during activity min ima as determined from the measured geo- 
magnetic aa-index (iLockwood et aI.L 12009 ). Upper panel: case 
without inflows. Lower panel: reference case with inflows. 
Numbers within symbols give the index of the cycle following 
the respective activity minimum. Full lines indicate regression 
curves, dashed lines the bisectrix of the first quarter plane. 



hemisphere, increasing the flux imbalance and thus eventually 
strengthening the axial dipole, and 2) the convergence of the in- 
flows reduces the tilt angle, counteracting the diffusive latitudi- 
nal spreading of a bipolar region, and thus weakening its even- 
tual contribution to the axial dipole. For weak to moderate cy- 
cles, the first effect dominates and the enhanced cross-equator 
transport leads to stronger axial dipole fields around activity 
minima. For stronger cycles, the combination of the Waldmeier 
effect (whereby stronger cycles peak earlier) and the fact that 
activity starts at higher latitudes before propagating towards the 
equator results in the activity belts being further away from the 
equator during the maximum phase of such cycles. This sup- 
presses the flux transport by cross-equatorial flows, so that even- 
tually the opposing effect 2) becomes dominant: the inflows act 
against the diffusion to reduce the tilt angle and thus reduce 
the diffusive cross-equatorial transport of magnetic flux. Fig. |2] 
shows this clearly for the minimum periods between cycles 19, 
20, and 21: while the maxima of the axial dipole flux in the case 
without inflows (red curve) largely follow the amplitudes of the 
preceding cycles in a linear manner, including the inflows (black 



3.2. Relation between axial dipole and cycle strength 

There is an empirical correlation between quantities related to 
the polar fields around activity minimum and the amplitude 
of the subsequent activity cycle. For instance, considering the 
values of the open heliospheric flux n ear activity minimum as 
infeiTed from geomagnetic variations (ILockwood et all l2009h 
and the maximum sunspot number (based on smoothed monthly 
mean valuefl) of the subsequent cycle, Rmax, we find r - 0.87 
(p = .0025) for cycles 13-21 . This is consistent with the re- 
sults of lWang & Sheelevl (l2009l) . On the other hand, the correla- 
tion with the maximum sunspot number of the preceding cycle 
is insignificant (r = 0.55, p = 0.12). These results indicate that 
the polar fields reflect the poloidal source of the toroidal field of 
the next cycle, which is consistent with the interpretation of the 
solar cycle in terms of a BL-type dynamo. 

In order to evaluate the effect of the latitudinal inflows on 
this coiTelation and thus infer their viability as a nonlinear feed- 
back mechanism for BL-type dynamos, we consider the coiTela- 
tions between the maxima of Odipoie as provided by SFT simu- 
lations and the values of R^nx of the subsequent and preceding 
cycles, respectively. Fig.|4]shows these relations for the simula- 
tions with and without inflows, respectively. In latter case, the 
coiTelation is moderate for Rmax of the preceding cycle (upper 
left panel; r = 0.68, p = 0.044) and insignificant for the sub- 
sequent cycle (upper right panel; r = 0.54, p = 0.13), which 
is in striking contrast to the empirical correlations. Including the 
inflows completely changes the picture: while the correlation of 
Odipoie with Rmax of the preceding cycle becomes insignificant 
(lower left panel; r = 0.38, p - 0.31), the correlation with R^ax 
of the subsequent cycle is highly significant in this case (lower 
right panel; r - 0.82, p = 0.0065), which is comparable to the 
empirical correlation based on the open heliospheric flux. 

2 .ftp : //ftp ■ ngdc . noaa . gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/doc; 
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Fig. 4. Correlation diagrams for the maxima of Odipoie during activity minima between cycles « - 1 and n from the SFT simulation 
vs. the maxima of the sunspot number of the the cycle preceding the minimum (n-1, left panels) and of the subsequent cycle (n, 
right panels), respectively. Upper row. case without inflows; lower row: reference case with inflows. Numbers within symbols give 
the index, n, of the cycle following the respective activity minimum. Linear regressions are indicated by the dashed lines. 



These results show that the inflows towards the activity belts 
provide an important nonlinearity in the evolution of the Sun's 
poloidal field at the solar surface. This conclusion follows essen- 
tially from the magnitude of the modeled inflows and is largely 
independent of the model parameters. Furthermore, our results 
strongly suggest that the latitudinal inflows towards the activity 
belts have a dominant influence on the strength of the subse- 
quent cycle and provide a nonlinear feedback mechanism for a 
BL-type dynamo. While other nonlinear processes might also 
play a role, the high correlations found here suggest that their 
effect is probably rather limited. Further work is going to assess 
the extent to which the good correlation with the observed cy- 
cle strengths found here depends on the SFT parameters and on 
our model of the inflows. Here we simply note that the param- 
eters u sed in this study are es sentially those of Cameron et al.l 
(l20Tol SFT parameters) and of lCameron & Schiissler! (1201 Oi in- 
flow model), which all are based on observational constraints. 



The high sensitivity of the axial dipoie strength on the am- 
pUtude of the inflows is consistent with the large variation of 
the cycle amplitudes in the historical record, indicating also that 
rather weak fluctuations of the surface distribution of active re- 
gions could temporarily switch off the BL-type dynamo and may 
drive the system into an extended minimum state. 



4. Conclusion 

Our SFT simulations show that magnetic-field-dependent latitu- 
dinal inflows converging towards the activity belts significantly 
affect the build-up of the polar field by modifying the cross- 
equator transport of magnetic flux. The resulting amplitudes of 
the magnetic flux contained in the axial dipoie component dur- 
ing activity minima correlate well with the empirically derived 
values of the open heliospheric flux during these periods. The 
inflows strengthen the axial dipoie in weaker cycles, for which 
parts of the inflows provide enhanced cross-equator transport of 
magnetic flux. For strong cycles, the reduction of the tilt angles 
of bipolar magnetic regions by the converging inflows dominates 
and leads to a weakening of the axial dipoie. Consistent with the 
empirical results, the SFT simulations including latitudinal in- 
flows show a strong correlation between the axial dipoie around 
activity minimum and the observed maximum sunspot number 
of the subsequent cycle. This indicates that the inflows are a key 
ingredient in determining the amplitude of solar cycles by pro- 
viding a nonlinear feedback mechanism for the saturation of a 
Babcock-Leighton-type dynamo mechanism. 
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